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Black hole candidates in X-ray binary systems and at the centers of galaxies are expected to be 
the Kerr black holes of General Relativity, but the actual nature of these objects has yet to be 
verified. In this paper, we consider the possibility that they are exotic compact objects and we 
describe their exterior gravitational field with a subclass of the Manko-Novikov metrics, which are 
exact solutions of the vacuum Einstein's equations and can describe the spacetime geometry around 
bodies with arbitrary mass-multipole moments. We point out that around a Manko-Novikov object 
there may exist many disconnected non-plunging regions at small radii, with no counterpart in 
the Kerr background, and that their existence may be tested. For instance, in the presence of an 
accretion disk, they may be filled by the accreting gas, forming a ring structure that might remind 
the rings of Saturn. We suggest that the existence of these regions may have a clear observational 
signature in the waveform of the gravitational radiation emitted by an EMRI: in the last stage of 
the inspiral, the waveform would be the combination of "regular chirps" , produced when the small 
object orbits in one of the non-plunging regions, and "bursts" , released when the small object jumps 
from a non-plunging region to another one at smaller radii. Our conclusions are supported by some 
numerical calculations of trajectories in the geodesic approximation, in which a particle plunges 
from the ISCO and then seems to get trapped in the potential well at smaller radii. 

PACS numbers: 97.60.Lf, 04.50.Kd, 97.10.Gz 



I. INTRODUCTION 

In 4-dimensional General Relativity, uncharged black 
holes (BHs) are described by the Kerr solution and are 
characterized by only two parameters: the mass, M , and 
the spin angular momentum, J [T]. Astrophysical BH 
candidates are stellar-mass compact objects in X-ray bi- 
nary systernsp and super-massive bodies at the centers 
of galaxies jSlM The former have a mass exceeding 3 Mq 
and are therefore too heavy to be neutron or quark stars 
for any plausible matter equation of state [S]. At least 
some of the super-massive BH candidates in galactic nu- 
clei turn out to be too heavy, compact, and old to be 
clusters of non- luminous bodies [B] . These two classes of 
objects are thought to be the Kerr BH predicted by Gen- 
eral Relativity because there is no alternative explanation 
in the framework of conventional physics. However, their 
actual nature has yet to be proven [7]. 

In this work, we consider the possibility that astrophys- 
ical BH candidates are exotic compact objects and we 
describe their exterior gravitational field with a subclass 
of the Manko-Novikov (MN) metrics [8j. The MN space- 
times are exact solutions of the vacuum Einstein's equa- 
tions and they are characterized by the mass of the com- 
pact object, its spin angular momentum, and an infinite 
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^ I'hc existence of a third class of objects, intermediate-mass BH 
candidates with AI ~ 10^ — 10* Mq, is still controversial, because 
their detections are indirect, and definitive dynamical measure- 
ments of their masses are lacking [4j . 



number of deformation parameters, setting all the mass- 
multipole moments. Previous work on these metrics con- 
sidered only the simplest case with one deformation pa- 
rameter, the anomalous mass-quadrupole moment q, as- 
suming zero all the higher anomalous mass moments [Sj- 
IT^ . The existence of some disconnected non-plunging 
regions at small radii in the MN spacetimes, which have 
no counterpart in the Kerr BH background (but see JJij ) , 
is already known. So far they have been quite ignored 
and were considered to be of no astrophysical interest. 
The reason is that it was thought that their orbital en- 
ergy is always higher than the one at larger radii and 
therefore that they would be difficult to populate by the 
gas of the accretion disk or by small bodies inspiralling 
into the BH candidate. Here we point out that this is not 
always true and that the existence of such non-plunging 
regions may have astrophysical implications. Moreover, 
in addition to a non-vanishing q, we consider the possi- 
bility of a non-vanishing anomalous mass-hexadecapole 
moment h, and we argue that the maximum number of 
non-plunging regions increases as the number of non-zero 
deformation parameters increases. 

In the Kerr metric, equatorial circular orbits are al- 
ways vertically stable, while they are radially stable for 
radii r > nsco, where nsco is the radius of the in- 
nermost stable circular orbit (ISCO). Equatorial circu- 
lar orbits with radii rpo < r < nscO; where rpo is 
the radius of the photon orbit, are unstable for small 
perturbations along the radial direction, while there are 
no equatorial circular orbits for r < rpo- In non-Kerr 
spacetimes, equatorial circular orbits may also be verti- 
cally unstable, which leads to a number of completely 
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new phenomena [HI [T3] . In the MN metric with a non- 
vanishing anomalous quadrupole moment, there may be 
an island of stable equatorial circular orbits in the region 
T < nsco; that is, there may exist two radii ri and r2 
such that equatorial circular orbits with ri < r < r2 are 
stable. If we consider the MN metric with non-vanishing 
anomalous mass-quadrupole and mass-hexadecapole mo- 
ment, there may exist two disconnected non-plunging re- 
gions for r < nscO: say ri < r < r2 and r4 < r < r^,. 
In the presence of an accretion disk, these regions may 
be filled by the accreting gas, because their energy and 
angular momentum can be lower than the ones at larger 
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Here R — + — 1 and P„ are the Legendre poly- 
nomials with argument xy/R, 

P - P 
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The MN solutions have an infinite number of free pa- 
rameters: a, which regulates the spin J of the central 
compact object; k, which regulates the mass M of the 
central compact object; and a„ {n = 1, -|-cx)) which 



radii. 

Throughout the paper, we use units in which Gn = 
c = 1. 



II. MANKO-NOVIKOV SPACETIMES 

The MN metrics [S] are stationary, axisymmetric, 
and asymptotically flat exact solutions of the vacuum 
Einstein's equations with arbitrary mass-multipole mo- 
ments. In quasi-cylindrical coordinates (p, z) and in pro- 
late spheroidal coordinates (x,?/), the line element is, re- 
spectively, 
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regulate the mass-multipole moments, starting from the 
dipole at, to the quadrupole a2, etc. For a ^ and 
a„ = 0, the metric reduces to the Kerr solution. For 
a = q;„ = 0, it reduces to the Schwarzschild one. For 
a = and a„ ^ 0, we obtain the static Weyl metric. 

Without loss of generality, we can put ai = to bring 
the massive object to the origin of the coordinate sys- 
tem. In what follows, we restrict our attention to the 
subclass of MN spacetimes with a„ = for n ^ 2 and 
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FIG. 1. The left panel shows the profile of the effective potential on the equatorial plane, Kff (p, z = 0), for M = 1, a* = 0.5, 
q = -0.18, h = -0.024, E = 0.9143, and = 2.83. Kff = (black horizontal line) is the limit of the allowed region. The 
roots of the effective potential, i.e., Aqb, Ajb, Cob, Cib, Eqb, and Ejb, determine the limits of three disconnected regions 
of bounded non-plunging orbits. The index OB stands for outer border, while IB for inner border. The minima j4niin, C'min, 
and -Emin of Vcff indicate the existence of stable circular orbits, while the maxima B and D indicate the existence of unstable 
circular orbits. The right panel shows the limits of the allowed region Vch = on the p ~ z plane, known as CZVs. 



^ We then have four free parameters {k, a, a-i, 04) 
related to the mass M, the dimensionless spin param- 
eter a.if = J/M"^, the anomalous mass-quadrupole mo- 
ment q ^ —{Q — QKcri)/M^, and the anomalous mass- 
hexadecapole moment h = — {H — Hq)/AI^, where Hq is 
the mass-hexadecapole moment of the object with h = 
and for q = it reduces to i?Kcrr — J'^/M^, by the rela- 
tions 
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Note that q measures the deviation from the quadrupole 
moment of a Kerr BH, but when g ^ also the higher- 
order mass-multipole moments have a different value 



^ A self-gravitating fluid is expected to be symmetric with respect 
to the equatorial plane. In this case, all the a„ with n odd must 
vanish and the n = 2 and 4 are the two possible independent 
leading order corrections to the Kerr geometry of the MN metrics 
(in the MN solutions, the current moments are not independent; 
e.g. the current-octupole moment is set by the mass-quadrupole 
moment). 



than in Kerr. For this reason, h measures the devia- 
tion from the hexadecapole moment of a Kerr BH only 
when g = 0. 

The relation between prolate spheroidal and quasi- 
cylindrical coordinates is given by 



p = k^/{x^ - 1) (1 - y^) , z^kxy, (15) 
with inverse 

x=^Up^ + {z + kf + ^p^ + {z- kf \ , 



y = ^[^Jp^ + iz + kf-^p' + {z-kf) ■ (16) 

The standard Boyer-Lindquist coordinates (r, 9) are re- 
lated to the prolate spheroidal coordinates and the quasi- 
cylindrical coordinates by 



COS0 = y , (17) 

p= - 2Mr + alM^ sin 6*, 

z = (r-M) COS 6*. (18) 



and 



Because the MN metrics are stationary and axisym- 
metric, geodesic orbits have two constants of motion, the 
specific energy E = —ut and the z-component of the 
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specific angular momentum = u^. The t- and 4>- 
components of tlie 4-velocity of a test-particle are there- 
fore 
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where the dots denote derivatives with respect to the 
proper time. From the normalization of the 4-velocity 
1, we can write 



p^ + P + V,s{E,L,p,z)=Q, 



where the effective potential V^ff is defined by 
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Due to stationarity and the axisymmetricity, the study 
of the geodesic dynamics can be restricted on a meridian 
plane {(j) = const), which is implied by Eq. (22) as well. 



Moreover, the motion takes place inside a region defined 
by Vcs < (left panel of Fig. [T]), whose border is the curve 
of zero velocity (CZV) Kg = (right panel of Fig.[l]). 

Circular orbits on the equatorial plane must have p — 
i = 0, which implies Vcs — 0, and p = z = 0, which im- 
plies dpVcs — and dzVcS — 0. This means that circular 
equatorial orbits are located at simultaneous zeros and 
extrema of the effective potential. The extrema shown in 
the left panel of Fig. [T] indicate the existence of these cir- 
cular orbits, the minima ^min? C*niim and -Emin indicate 
the existence of stable circular orbits, while the maxima 
B and D indicate the existence of unstable circular or- 
bits. In order to calculate the positions of the circular 
orbits, we must also demand that the effective potential 
to be equal to zero. Thus, because dzVeS = is satis- 
fied identically for z = (simply because of the reflection 
symmetry of our MN metric with respect to the equato- 
rial plane), from the aforementioned conditions one can 
obtain E and as a function of the radius of the orbit 
P- 



9tt + gt4,^ 
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is the angular frequency and the sign -I- (— ) is for or- 
bits with angular momentum parallel (antiparallel) to the 
spin of the compact object. These orbits are stable under 
small perturbations in the radial direction if dpVcs > 0, 
and in the vertical direction if d^Vcs > 0. 
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FIG. 2. The region with CTCs (in gray) is determined by 
the condition g^cj, < 0, while the regions where motion is 
permitted are delimited by the CZV (black curves). It seems 
that the two regions never overlap for plausible values of the 
orbital parameters. The MN parameters, E, and Lz are the 
same as in Fig. [l] 



The no-hair theorem jTj states that the only asymp- 
totically flat, vacuum and stationary solution of the Ein- 
stein's equations that is non-singular on and outside an 
event horizon and that presents no closed time- like curves 
(CTCs) outside it is given by the Kerr metric. Therefore, 
the MN spacetimes must either have no event horizon, or 
present naked singularities or CTCs outside it. In fact, a 
MN metric can present the last two anomalies. Namely, 
the surface x = 1 {p = 0, \z\ < k), which is the event 
horizon in the Kerr case a„ = 0, is in general (a„ 7^ 0) 
only a partial horizon, because it presents a naked cur- 
vature singularity on the equatorial plane (i.e. at a; = 1, 
y = 0, corresponding to p = 2: = 0). Moreover, there 
are CTCs outside it (Fig. [2|. However, it seems that the 
CTCs and the allowed region Vctf < do not overlap. 

Nevertheless, in our scenario the exotic object is com- 
pact enough so that the plunging particles can reach the 
inner stable orbits, but not so much to have the region 
with CTCs. Furthermore, the pathological region with 
CTCs can even be removed by exciting the spacetime. 
This occurs, for example, if the spacetime goes to a new 
phase and a domain wall is formed [HI [16] . Across the 
domain wall, the metric is nondifferentiable and the ex- 
pected region with closed timelike curves arises from the 
naive continuation of the metric ignoring the domain 
wall. The latter can be made of very exotic stuff, and 
in the literature there are examples of supertubes |15| 
and fundamental strings |16| . 
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III. PROPERTIES OF EQUATORIAL 
CIRCULAR ORBITS IN MN SPACETIMES 



A. MN spacetimes with g 7^ and h = 

This is the case aheady discussed in the hterature. The 
basic properties of equatorial circular orbits can be un- 
derstood by considering a specific spin parameter a* and 
by changing the anomalous quadrupole moment q. As an 
example, we report the values for the case — 0.5 (see 
Tab.|l]). 

For q = 0, we recover the Kerr spacetime. Outside the 
even horizon, the spacetime is everywhere regular (there 
are no CTCs). Circular orbits on the equatorial plane 
exist for radii p > ppo, they are always vertically stable, 
but they are radially stable only for p > pisco- If the 
compact object is surrounded by a thin accretion disk, 
the inner edge of the disk is at pisco- For a* — 0.5, 
we have ppo — 1.03 and pisco — 3.11 (here and in the 
following, we use units in which Af — 1). 

MN spacetimes with q ^ have regions with CTCs 
just outside the event horizon. That is true even on 
the equatorial plane. When the compact object is more 
oblate than a Kerr BH (g > 0), there are no interesting 
new properties with respect to the Kerr metric. As the 
object becomes more and more oblate {q increases), the 
value of Pisco increases. For the case q ~ 0.20 shown in 
Tab. |T] equatorial circular orbits exist for p > 1.39, they 
are always vertically stable, but they are radially stable 
only for p > pisco — 3.50. Unlike the Kerr background, 
a narrow region with circular orbits appears even at very 
small radii (p ~ 0.04 for q — 0.20). This narrow region 
extends up to the region with CTCs. Even if it seems 
that this region never overlaps with the CTC region, it is 
more likely unphysical: these orbits have negative energy, 
suggesting that their possible existence would make the 
object unstable [T7]. While orbits with negative energy 
exist even in the ergoregion of a Kerr BH, they are never 
stable. 

For g < 0, the spacetime can have properties absent 
in the Kerr background and that are potentially of as- 
trophysical interest. Let us consider again the exam- 
ple with a* — 0.5. For small \q\, there is a region with 
circular orbits at small radii. Such orbits are radially 
stable, but vertically unstable. If we consider the case 
q = -0.01 of Tab. [T) such a region is at 0.25 < p < 0.35. 
At larger radii, circular orbits exist for p > 0.99, they 
are always vertically stable, but they are radially stable 
only for p > pisco = 3.09. As |g| increases, the size 
of the region with circular orbits at small radii increases 
as well, and eventually it merges with the one at larger 
radii. At this point, we have two disconnected regions 
with stable orbits, Pi < p < P2 and p > pisco! this is 
the case oi q — —0.1, with pi = 1.00, p2 — 1.16 and 
PISCO = 2.83, and of q -0.2, with pi = 1.24, pa = 2.00 
and Pisco = 2.17. The orbit with radius p = pi is 
marginally unstable along the vertical direction, while 
the one with p = p2 is marginally unstable along the ra- 
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FIG. 3. MN spacetimes with non-vanishing anomalous 
quadrupole moment q. The solid red curves separate space- 
times with qualitatively diflerent properties. Region 1: space- 
times in which equatorial circular orbits of astrophysical in- 
terest are located at p > pisco and the orbit at p = pisco is 
marginally stable along the radial direction. Region 2: space- 
times in which equatorial circular orbits of astrophysical inter- 
est are located at p > pisco and at pi < p < ps (ps < pisco); 
the orbit at p = pi is marginally stable along the vertical 
direction, while the one at p — pisco is marginally stable 
along the radial direction; the energy of the orbits at p = ps 
and a,t p — pisco is the same and therefore the gas's parti- 
cles can plunge from the inner edge of the external disk to 
the outer edge of the internal disk. Region 3: spacetimes in 
which equatorial circular orbits of astrophysical interest are 
located at p > pi and the orbit aX p = p\ is marginally sta- 
ble along the vertical direction. The compact object is oblate 
if g > — aj; the green dashed curve separates oblate objects 
from prolate objects. The compact object is more oblate than 
a Kerr black hole if g > 0. See the text for details. 



dial direction. If the orbits in the region pi < p < P2 
have energy larger than the orbit at p — pisco (the 
case q = —0.1), the region pi < p < p2 is not of as- 
trophysical interest, as particles/gas orbiting around the 
compact object after reaching the radius pisco plunge 
onto the compact object, without populating the region 
pi < p < p2- However, as \q\ increases, the size of the 
region pi < p < P2 also increases and there may exist 
a region pi < p < p^ with pa < p2 in which the or- 
bital energy is smaller than the one at p = Pisco (the 
case q = —0.2). These orbits are of astrophysical inter- 
est because they can be filled up by the gas of accretion 
reaching the radius p ~ pisco- If 1 9 1 continues increas- 
ing, the region P2 < p < Pisco with radially unstable 
orbits disappears: now circular orbits are always radially 
stable, but they are vertically stable for p > pi, as is for 
the case q = —0.3, with pi = 1.43. 

Fig. [3] shows the three different kinds of MN space- 
times with q 7^ on the plane spin parameter-anomalous 
quadrupole moment. The region 1, which includes the 
Kerr solution q = 0, has equatorial stable circular orbits 
for p > Pisco and the orbit at p = pisco is marginally 
radially stable; for q < 0, these spacetimes may also 
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have equatorial stable circular orbits at smaller radii, but 
they are not of astrophysical interest because these or- 
bits require energies higher than the one of the orbit at 
P = Pisco- In the region 2, there are spacetimes with 
equatorial stable circular orbits of astrophysical interest 
for pi < p < P3 and p > pisco {ps < Pisco); the orbit 
at p = pi is marginally vertically stable, while the one at 
P = Pisco is marginally radially stable. Lastly, the region 
3 includes those spacetimes in which equatorial circular 
orbits are always radially stable, but they are vertically 
stable only for p > pi. 



B. MN spacetimes with q — and h 

If the deformation parameter of the spacetime is the 
anomalous hexadecapole moment h instead of the q, the 
picture is qualitatively the same as the one discussed in 
the previous subsection, with h playing the role of — g 
(the minus sign is just a matter of definition). 



C. MN spacetimes with q < and h 

Let us now discuss the more general scenario, when 
both q and h may be non-vanishing, starting from the 
case with g < 0. As in Subsection [III A[ we can consider 
the specific case a* = 0.5 (see Tab. |II[); for another value 
of a*, the picture is qualitatively the same, but the same 
phenomena occur for different values of q and h. First, 
we consider a spacetime in which \q\ is large enough that 
for /i = we have stable equatorial orbits for p > pi. In 
Tab. [nj we consider the specific case q = —0.25, where 
Pi — 1.34 when h = 0. If we change the anomalous hex- 
adecapole moment to h = —0.03, we find that circular 
orbits exist even in two narrow region at small radii: a re- 
gion close to the one with CTCs (p « 0.27), and another 
one in which the orbits are vertically stable but radially 
unstable (0.40 < p < 0.43). The latter expands as \h\ 
increases and it merges with the region of circular orbits 
at large radii. For instance, the regions 0.40 < p < 0.43 
and p > 0.68 for h = —0.03 become the region p > 0.53 
when h = —0.05. At this point, we have three regions of 
stable circular orbits: one close to the pathological space 
with CTCs, one at p4 < p < ps, and the last one for 
p > pi. The orbit at p = p^ is marginally radially stable, 
while the one at p — p^ is marginally vertically stable. 
For h = —0.05, the region close to the one with CTCs is 
at p w 0.32, while p4 = 0.67, pg = 0.70, and pi = 1.17. 
A spacetime with similar properties does not exist if the 
only non-vanishing deformation parameter is q. Let us 
notice that the region p^ < p < p^ may include a subre- 
gion of astrophysical interest p4 < p < p^, where pe < Ps, 
in which the particles have energies smaller than the one 
of a particle at p = pi. In presence of an accretion disk, 
the region p4 < p < pg could be filled by the accreting 
gas. For larger values of \h\, the gap ps < p < pi disap- 
pears and stable orbits exist for p > p^. This is the case 



with h = —0.08, where p4 = 0.91. By continuing to in- 
crease the value of a new gap shows up, in which the 
orbits are vertically stable but radially unstable. So, for 
h = —0.1 stable orbits exist for p^ < p < pj and p > ps, 
where p4 = 1.14, py — 1.71, and ps — 2.04. The new gap 
is p7 < p < p8 and we have used a different notation with 
respect to the gap ps < p < pi of the case h = —0.05 be- 
cause it has different properties: py and pg are marginally 
radially stable, while ps and pi were marginally vertically 
stable. If we continue increasing the value of \h\, the gap 
Pi < p < Ps disappears. This is the case h — —0.15: 
the spacetime now has features similar to that of a MN 
background with q > and vanishing h, in which stable 
orbits exist for p > p4 (p4 = 2.32 for h — —0.15), the 
orbit at the radius p = p4 is marginally radially stable, 
and the value of p4 increases as \h\ increases. 

If we start with q = —0.25 and h = 0.0 and we increase 
the value of h (now ft. > 0), we simply move pi at larger 
radii. For example, for ft = 1.0 we have pi = 2.06, while 
for ft — 5.0 we have pi = 2.83. We never form islands of 
stable circular orbits at small radii. 

The properties of the MN spacetime may be even more 
interesting when q is still negative, but |g| is not so large, 
so that for ft = the region pi < p < p2 is still separated 
by a gap from the one p > pisco- As a specific example, 
we have q = —0.15 in Tab. [llj The properties of the 
equatorial circular orbits can be understood on the basis 
of the previous discussion. When ft = —0.01, pi = 1.08, 
P2 = 1.40, and pisco — 2.64. For smaller ft, we may have 
three disconnected regions of stable circular orbits (and, 
at least in some cases, of astrophysical interest): a new 
region at pg < p < pig, the regions Pi < p < P2, and 
the one at p > pisco- In the specific case ft = —0.029, 
we find p9 = 0.68, pio = 0.77, pi = 0.89, p2 = 1.22, and 
Pisco = 2.67. 



D. MN spacetimes with q > and h 

The properties of the case q > can be easily expected 
on the basis of the results found in the previous subsec- 
tions. The effect of increasing/decreasing ft is equivalent 
to that of decreasing/increasing q. Tab. Ill shows the 
case q = 0.25. 



E. Islands of stable equatorial circular orbits of 
astrophysical interest 

The study of the previous subsections suggests the fol- 
lowing conjecture: the MN spacetimes with arbitrary de- 
formation parameters {q, ft, and higher order deforma- 
tions) can potentially have several disconnected regions 
with stable equatorial orbits of astrophysical interest. 
Accreting compact objects described by the MN solu- 
tions, if they exist, might thus be surrounded by Saturn- 
like structures, in which the accreting gas can fill different 
rings around the compact object. In the case a* = 0.5, 
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a* 


g 


h 


Circular orbits 


9pVeff > 


9?14ff > 


CTCs 


0.5 








p > 1.03 


p > 3.11 


p > 1.03 


no CTCs 


0.5 


0.20 





p w 0.04 


p « 0.04 




p < 0.04 








p > 1.39 


p > 3.50 


p > 1.39 




0.5 


-0.01 





0.25 < p < 0.35 


0.25 < p < 0.35 




p < 0.15 








p > 0.99 


p > 3.09 


p > 0.99 




0.5 


-0.10 





p > 0.56 


0.56 < p < 1.16, p > 2.83 


p > 1.00 


p < 0.35 


0.5 


-0.20 





p > 0.73 


0.73 <p< 2.00, p > 2.17 


p > 1.24 


p < 0.45 


0.5 


-0.30 





p > 0.85 


p > 0.85 


p > 1.43 


p < 0.52 



TABLE I. Regions with equatorial circular orbits (fourth column), radially stable orbits (fifth column), vertically stable orbits 
(sixth column), and CTCs (seventh column) in MN spacetimes with spin parameter a* = 0.5, anomalous mass-quadrupole 
moment q, and vanishing anomalous mass-hexadecapole moment h. See the text for details. 
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h 


Circular orbits 


d'^pVeH > 


9f V^ff > 


CTCs 


0.5 


-0.25 


0.00 


p > 0.79 


p > 0.79 


p > 1.34 


p < 0.49 


0.5 


-0.25 


-0.03 


p«0.27 
0.40 < p < 0.43 
p > 0.68 


p w 0.27 
p > 0.68 


p« 0.27 
0.40 < p < 0.43 
p > 1.26 


p < 0.27 


0.5 


-0.25 


-0.05 


p « 0.32 
p > 0.53 


p« 0.32 
p > 0.67 


ps:; 0.32 
0.53 < p < 0.70, p > 1.17 


p < 0.32 


0.5 


-0.25 


-0.08 


p 0.38 
p > 0.69 


p « 0.38 
p > 0.91 


p > 0.69 


p < 0.38 


0.5 


-0.25 


-0.10 


p« 0.41 
p > 0.77 


p«0.41 
1.14 < p < 1.71, p > 2.04 


p« 0.41 
p > 0.77 


p < 0.41 


0.5 


-0.25 


-0.15 


p w 0.45 
p > 0.92 


p w 0.45 
p > 2.32 


p > 0.92 


p < 0.45 


0.5 


-0.25 


1.00 


p > 1.35 


p > 1.35 


p > 2.06 


p < 0.87 


0.5 


-0.25 


5.00 


p > 1.88 


p > 1.88 


p > 2.83 


p < 1.19 


0.5 


-0.15 


-0.010 


p w 0.21 
0.30 < p < 0.31 

p > 0.59 


p w 0.21 
0.59 < p < 1.40, p > 2.64 


0.30 < p < 0.31 
p > 1.08 


p < 0.21 


0.5 


-0.15 


-0.029 


p « 0.30 
p > 0.54 


p « 0.30 
0.68 < p < 1.22, p > 2.67 


0.54 < p < 0.77, p > 0.89 


p < 0.30 


0.5 


-0.15 


-0.035 


p « 0.32 
p > 0.61 


p « 0.32 
0.80 < p < 1.13, p > 2.67 


p > 0.61 


p < 0.32 


0.5 


-0.15 


-0.050 


p w 0.36 
p > 0.74 


p w 0.36 
p > 2.69 


p > 0.74 


p < 0.36 



TABLE U. Regions with equatorial circular orbits (fourth column), radially stable orbits (fifth column), vertically stable orbits 
(sixth column), and CTCs (seventh column) in MN spacetimes with spin parameter a* = 0.5, anomalous mass-quadrupole 
moment q < 0, and anomalous mass-hexadecapole moment h. See the text for details. 



q = —0.15 and h = —0.029, we can have an accretion 
disk with two internal rings. However, for non-vanishing 
higher order deformations, three or more rings may be 
possible. 



IV. OBSERVATIONAL SIGNATURE 

Let us consider a spacetime in which equatorial circular 
orbits are stable for p > pisco in the inner regions 
Pi < p < p2, P4 < p < P5 , etc. at smaller radii. In the 
presence of a gas of accretion, we may have the usual thin 
disk at p > Pisco and several gas rings for pi < p < p^, 
P4 < p < Pq, etc. Is the thermal spectrum of this multi- 
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g 


h 


Circular orbits 


a^Kft > 


9f Kft > 


CTCs 


0.5 


0.25 


0.0 


p ^ 0.42 


0.42 




p < 0.42 








p > 1.44 


p > 3.58 


p > 1.44 




0.5 


0.25 


-1.0 


0.77 


0.77 


p^ 0.77 


p < 0.77 








p > 1.75 


p > 3.79 


p > 1.75 




0.5 


0.25 


0.1 


0.59 < p< 0.67 


0.59 < p< 0.67 




p < 0.42 








p > 1.37 


p > 3.55 


p > 1.37 




0.5 


0.25 


0.5 


p > 0.96 


0.96 < p< 1.55, p > 3.43 


p > 1.43 


p < 0.67 


0.5 


0.25 


1.0 


p > 1.17 


1.17 < p< 2.13, p > 3.21 


p > 1.73 


p < 0.79 


0.5 


0.25 


2.0 


p > 1.40 


p > 1.40 


p > 2.07 


p < 0.93 



TABLE III. Regions with equatorial circular orbits (fourth column), radially stable orbits (fifth column), vertically stable orbits 
(sixth column), and CTCs (seventh column) in MN spacetimes with spin parameter a, — 0.5, anomalous mass-quadrupole 
moment q > 0, and anomalous mass-hexadecapole moment h. See the text for details. 



part accretion disk very different from the one we can 
expect around a Kerr BH with inner edge Pjgco'^ The 
answer is no and the reason is the following. The high- 
energy part of the thermal spectrum of a thin accretion 
disk around a compact object depends on the radiative 
efficiency r; = 1 — E-m, where E-m is the specific energy 
of the gas particles at the inner edge of the disk. So, it 
does not matter if the disk has some gaps. The thermal 
spectrum of our disk will look like the one around a Kerr 
BH with the same radiative efficiency. As the gas in inner 
rings necessarily has energy lower than the one at pisco j 
the radiative efficiency of these objects is higher than 
r] = 1 — £'isco I which might be (erroneously) interpreted 
as a higher value of the spin parameter of the object if 
we do not take the radiation from the inner disk into 
account. The ring structure of the accretion disk may be 
observed in high-resolution images of the disk itself, see 
e.g. Figs. 5 and 6 of Ref. [TH], but this is definitively out 
of reach for present and near future facilities. 

How can we then identify the existence of islands of 
stable orbits around a BH candidate? We think that 
a clear observational signature of these non-Kerr space- 
times may be found in the gravitational wave signal emit- 
ted by an extreme-mass ratio inspiral (EMRI) system, 
in which a stellar-mass compact object orbits around a 
super-massive BH candidate. In the standard case, the 
frequency and the amplitude of the waveform of the grav- 
itational radiation are expected to increase regularly as 
the orbit of the small body shrinks to the ISCO radius 
- this is usually called the "chirp" phase in the litera- 
ture. At the ISCO, the orbit becomes radially unstable 
and the stellar-mass body plunges into the super-massive 
BH candidate within a dynamical time-scale, emitting 
an irregular waveform. In the presence of inner regions 
of stable orbits, the picture may change as follows. As 
in the Kerr metric, the frequency and the amplitude of 
the waveform of the gravitational radiation should in- 
crease regularly as the small body approaches the ISCO 
radius. However, once at the ISCO, the small body may 
not plunge into the central object, but rather land in 



an inner region with stable orbits. During the transi- 
tion from the ISCO to the island, the system would emit 
some complicated waveform of gravitational waves to ad- 
just the small body to the new orbit, but, after that, we 
should start observing again the typical EMRI waveform, 
in which the frequency and the amplitude increase regu- 
larly. This process may be repeated several times, with 
the small body that reaches the inner edge of an outer 
region and jumps to another region at smaller radii. The 
whole waveform should thus consist of several pieces of 
regular waveforms (the usual EMRI-type signal) broken 
by "bursts" , due to gravitational waves emitted when the 
small body adjusts into the new orbit. 

Of course, our scenario can effectively work if the small 
body does jump from one region to another, without 
plunging directly into the central object. The correct 
approach to see if this is indeed the case is to compute 
the trajectory of one of these bodies including the radi- 
ation reaction. However, there is not an established pre- 
scription in the bibliography (at least to our knowledge) 
for handling radiation reactions in non-Kerr spacetimes. 
There are some tricks one can use in the relatively weak 
field limit [H [10] , but they cannot be used here, as in 
our case we are considering the region close to the com- 
pact object, where gravity is strong. Moreover, when the 
small body leaves the inner edge of a region of stable or- 
bits, its trajectory is surely neither circular nor eccentric, 
and there might even be a chance that the orbit gets non- 
equatorial due to vertical instabilities, which makes the 
issue of radiation reaction even more difficult to handle in 
the strong regime. Therefore, in support of our scenario, 
we have simply checked if (in the geodesic approximation) 
the small body plunges quickly into the central one af- 
ter having reached the ISCO radius or if it gets trapped 
in an inner region formed around circular orbits. Our 
numerical calculations show that the small body would 
indeed be trapped in an inner region, without plunging 
directly to the central object. 

In the specific example we examined, we used a MN 
spacetime with parameters M = 1, a* = 0.5, q = —0.15, 
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FIG. 4. The left panel shows the effective potential along p on the equatorial plane {z = 0) for the MN spacetime with 
parameters M = 1, a, = 0.5, q = -0.15, and h = -0.029. The orbital parameters are E = 0.911388 and = 2.820706, 
corresponding to the ones of a circular orbit at the ISCO radius. For the same spacetime and orbital parameters, in the right 
panel we show the CZVs (black curves) and an equatorial orbit {z = z = 0) leaving the ISCO (the gray curve) . In the embedded 
plot, the gray curve shows a non-equatorial orbit leaving the ISCO. 



h = —0.029. The effective potential on the equatorial 
plane for energy and angular momentum corresponding 
to a circular orbit at the ISCO radius (left panel of Fig. Ill) 
shows that between the ISCO and the horizon {p = 0} 
lies a deep potential well, which indicates the existence 
of stable circular orbits. On the p — z plane in the right 
panel of Fig.|4j the CZVs (black curves) show the allowed 
region for an orbit leaving the ISCO. The deep well of the 
Kff corresponds to the loop lying in the 0.3 < p < 0.6 
interval. If we evolve an equatorial geodesic orbit starting 
near p = pisco , the orbit goes straight to the center of the 
well (gray curve in the right panel of Fig.[4|, which means 
that it will be trapped for a while in the inner region of 
circular orbits, before the final plunge. If we now evolve 
a non-equatorial geodesic orbit starting near the ISCO 
radius, the orbit oscillates around the equatorial plane, 
but it also ends up in the well (embedded plot in the right 
panel of Fig. |4|. Thus, even some vertical instabilities 
may not cause a direct plunge and may allow the small 
body to jump to an inner region of stable circular orbits. 
In the Kerr space-time, the particle cannot get trapped 
at small radii simply because there is no deep potential 
well between the object and the ISCO. 

If we reduce the levels of energy and angular mo- 
mentum, as would happen due to the radiation reac- 
tion, we find two marginally stable circular orbits, see 
Fig. [5] The first one is radially stable and vertically 
unstable and has orbital parameters E = 0.888672 and 
Lz — 2.773112. The second one is radially unstable and 
vertically stable, with orbital parameters E = 0.796835 



and Lz = 2.378980. As the stellar-mass body leaves 
the ISCO radius, it loses energy and angular momentum. 
This changes the CZVs. In particular, the large allowed 
region in Fig. |4] may split. Depending on how the small 
body loses energy and angular momentum, it may adjust 
in the inner stable region at larger or smaller radii. 



V. SUMMARY AND CONCLUSIONS 

Astrophysical BH candidates are thought to be the 
Kerr BHs predicted by General Relativity, but this con- 
clusion is not yet supported by robust observational ev- 
idence. In this paper, we discussed the observational 
implications in the case that BH candidates are exotic 
compact objects whose exterior gravitational field is de- 
scribed by the Manko-Novikov spacetimes, which are ex- 
act solutions of the vacuum Einstein's equations. Such a 
possibility is not in contradiction with current observa- 
tions, but it requires some form of exotic matter beyond 
standard physics. However, the basic scope of our work 
is to provide an observational test of whether the centers 
of galaxies are occupied by Kerr BHs or by something 
more exotic, and not to argue that these exotic objects 
exist. 

Equatorial circular orbits around a Kerr BH are always 
vertically stable, while they are radially stable only for 
radii larger than nsco- The existence of stable equato- 
rial circular orbits at radii smaller than rigco in non-Kerr 
spacetimes was already known, but so far they have never 
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FIG. 5. The CZVs for the same spacetime parameters of Fig. [4]but for different orbital parameters. The left panel corresponds 
to = 0.888672 and = 2.773112, and the right to E = 0.796835 and = 2.378980. 



been considered to be of astrophysical interest as it was 
thought that their orbital energy was higher than the one 
at the ISCO and therefore that they could not be filled 
by the accretion gas or by small bodies inspiralling into 
the compact object. In this paper, we have pointed out 
that this is not always true. We have also argued that 
the maximum number of these non-plunging regions may 
increase as the structure of the spacetime becomes more 
and more complex by adding higher anomalous mass mo- 
ments. Lastly, we have proposed that the existence of 
these regions of stable equatorial circular orbits around 
astrophysical BH candidates may be tested by observing 
the gravitational waves emitted by a stellar-mass com- 
pact object inspiralling into a super-massive BH candi- 
date. In the standard case of a Kerr BH, the frequency 
and the amplitude of the waveform should regularly in- 
crease, in the so-called "chirp" of the system, up to the 



plunging of the small body into the BH. In the presence of 
these inner regions, the total waveform could instead con- 
sist of "regular chirps", produced when the stellar-mass 
compact object orbits in one of these stable regions, and 
"bursts" , necessary to adjust the small body to a new sta- 
ble orbit of a more internal region. Such a signature in 
the waveform might be observed with future space-based 
gravitational wave detectors. 
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